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Abstract 

By using a reductive perturbation technique applied to a two-level model, a generic two- 
dimensional sine-Gordon evolution equation governing the propagation of femtosecond spatiotem- 
poral optical solitons in Kerr media beyond the slowly-varying envelope approximation is put 
forward. Direct numerical simulations show that, in contrast to the long- wave approximation, no 
collapse occurs, and that robust (2+l)-dimensional ultrashort light bullets may form from ade- 
quately chosen few-cycle input spatiotemporal waveforms. In contrast to the case of quadratic 
nonlinearity, the light bullets oscillate in both space and time, and are therefore not steady-state 
lumps. 

PACS numbers: 42.65.Tg, 42.65.Re, 05.45.Yv 
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I. INTRODUCTION 



Since the first experimental realization more than one decade ago of two-cycle optical 
pulses in mode-locked Ti-sapphire lasers using double-chirped mirrors the field of few- 

cycle pulses (FCPs) has grown into one of the major area of modern ultrafast optics. These 
studies led to remarkable achievements in the field of extreme nonlinear optics, time-resolved 
laser spectroscopy, generation of soft X-ray radiation and isolated and controlled attosecond 
light pulses in the extreme ultraviolet regime, thereby opening the door for real-time probing 
of fast electron dynamics associated with excited-state atomic and molecular processes 4|- 
10| . These FCPs are key elements to many time-domain applications, driving the quest for 
stable mode-locked laser sources with ultrahigh spectral bandwitdth. A remarkable recent 
advance in this area was the syntesis of a sing le cycle of light in the near-infrared regime with 
compact erbium-doped fiber technology ll|. Another recent achievement is the possibility 
to use different optical parameter amplifiers to produce FCPs over a broad frequency range. 
The optical parametric chirped pulse amplification method has become one of the standard 
procedures to generate high peak power pulses with durations of a few optical cycles. Two- 



and three-cycle pulses with carrier wave 
mid-infrared were recently generated [12 



engths almost continuously tunable from visible to 



13|. 



On the theoretical arena the studies of t 
main directions: (i) the quantum approach 



le unique physics of FCPs concentrated on three 
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ITj , (ii) the refinements within the framework 



of the slowly varying envelope approximation (SV^ 



velope equations 
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2l|, and non-SVEA models 



EA 


of 


22 


26 



the nonlinear Schrodinger-type en- 



2q |. In Kerr-like media the physics of 



(l-l-l)-dimensional FCPs can be adequately described beyond the SVEA by using different 
dynamical models, such as the modified Korteweg-de Vries (mKdV) 2^, sine-Gordon (sG) 
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24| , or mKdV-sG equations 25|, |26| . The above mentioned evolution equations admit 



breather solutions, which are most suitable for describing the physics of few-optical-cycle 
solitons. Notice that other non-SVEA models 27H29|. especially the so-called short-pulse 
equation 30||, have been proposed. Moreover, a non-integrable generalized Kadomtsev- 
Petviashvili (KP) equation 3l| (a two-dimensional version of the mKdV model) was also 



put forward for describing the (2+l)-dimensional few-optical 



propagation in cubic nonlinear media beyond the SVEA 
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^cle spatiotemporal soliton 



33| . However, it has recently 



been shown that collapse of FCPs occurs for high enough intensities in the frame of cubic 



2 



generalized Kadomtsev-Petviashvili model, which describes ultrashort spatiotemporal op 



tical pulse propagation in cubic (Kerr-like) media 3J]. Notice that of particular interest 
for the physics of the (l+l)-dimensional FCPs is the mKdV-sG equation; thus by using a 



25| that the propagation of ultrashort 



system of two-level atoms, it has been shown in Ref. 
pulses in Kerr optical media is fairly well described by a generic mKdV-sG equation, which 
was also derived and studied in Refs. 



35 



36|. Recently, we showed that the mKdV-sG 



model is the most general of all approximate non-SVEA models for FCPs, and in fact con- 



tains all of them 



37|. Also, the propagation of single-cycle gap solitons in subwavelength 



periodic structures, without the use of the SV. 
been theoretically predicted in a recent study 



A of the Maxwell-Bloch equations has also 
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As concerning the physical assumptions made in order to get the above mentioned generic 
nonlinear evolution equations we first mention that soliton propagation implies that damping 
can be neglected. In dielectric media, this occurs far from any resonance frequency. Let us 
consider a two-level model with characteristic frequency Q, and denote by a frequency 
characteristic for the FCP soliton under consideration. The transparency condition implies 
that either Uw <^ Q or Q Uu,. The former case ^ ^) corresponds to the long wave 
approximation, whereas the latter case corresponds to the short wave approximation. In the 
frame of a two-level model, a mKdV equation is obtained if the frequency of the transition 
Q is far above the characteristic wave frequency (the long- wave approximation regime), 
while a SG model is valid if Q is much smaller than (the short-wave approximation 
regime). However, some of the transitions frequencies Qj are much smaller than u^, and 
the other ones much larger than Uw. This more realistic situation can be modeled by 
considering two transitions only, with different frequencies Qi and 1^2. The physical system 
is thus equivalent to a two-component medium, each of the two components being described 
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26 



35 



by a two-level model. As a result, a mKdV-sG model was put forward 
which is completely integrable in certain particular cases by means of the inverse scattering 
transform. It admits stable solutions of 'breather' type, which also give a good description 
of few-optical-cycle soliton propagation. 

The propagation of FCPs in a quadratic medium has also been described by either a 
KdV or a KP equation, in (1+1) or (2+1) dimensions, respectively, which evidenced either 
the stability of a plane wavefront, for a normal dispersion, or the formation of a localized 
spatiotemporal half-cycle soliton, for an anomalous dispersion 
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40| . By using a multi- 
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scale analysis, a generic KP evolution equation governing the propagation of femtosecond 
spatiotemporal optical solitons in quadratic nonlinear media beyond the SVEA was recently 
put forward |40[. Direct numerical simulations showed the formation, from adequately cho- 
sen few-cycle input pulses, of both stable line solitons (in the case of a quadratic medium 
with normal dispersion) and of stable lumps (for a quadratic medium with anomalous dis- 
persion). The perturbed unstable line solitons decay into stable lumps for a quadratic 



nonlinear medium with anomalous dispersion j40|. However, in Ref. 40|] we considered a 
set of two-level atoms and we have assumed that the characteristic frequency ujw of the con- 
sidered electromagnetic wave in the optical spectral range is much less than the transition 
frequency Vt of the atoms, i.e., we worked in the so-called long- wave approximation regime. 

Summarizing, in the long wave approximation, for a quadratic nonlinearity, half-cycle 
light bullets in the form of a single lump may exist, while for a cubic nonlinearity collapse 
occurs. In the present work we will show that in the short wave approximation for a cubic 
nonlinearity, a third physical situation occurs: few-cycle light bullets may form, oscillating 
in both space and time. 

More precisely, the aim of this paper is to derive a generic partial differential equation 
describing the dynamics of (2+l)-dimensional spatiotemporal solitons beyond the SVEA 
model equations in the so-called short-wave approximation regime, when we consider the 
situation in which the resonance frequency Vt of the two-level atoms is below the character- 
istic optical frequency Uyj. Several generalizations to (2+1) dimensions of the sG equation 
exist in the literature and can be derived in the short-wave approximation; some of them 
support localized solitons which are not oscillating as the sG breather and one-dimensional 
FCP solitons in Kerr media 
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44j |. The question whether the (2+l)-dimensional general- 



ization of sG equation, which is valid for FCP propagation in cubic nonlinear media, is one 



of the nonlinear dynamical systems discussed in Refs. |4l|-|44| or some other one, can only 
be decided by means of a rigorous derivation startin g fro m the basic equations. We use a 
reductive perturbation method (multiscale analysis) 45| and as a result we get a generic 
two-dimensional sine-Gordon (2D sG) equation to describe the propagation of ultrashort 
(2+l)-dimensional light bullets in a system of two- level atoms beyond the SVEA. This is 
the first result of the present work. 

This work is organized as follows. In Sec. II we derive the 2D sG equation governing 
the propagation of femtosecond spatiotemporal solitons beyond the SVEA. Direct numerical 
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simulations of the 2D sG nonlinear partial differential equation are given in Sec. Ill, where 
it is shown the generation of stable (2+l)-dimensional FCPs from perturbed FCP plane 
waves. This is the second result of this paper. We note that a numerical resolution of 



the 2D sG equation was published in an earlier work [46l], which showed the existence of 
ultrashort light bullets. However, it was believed that their amplitudes were not constant 
but decreased during propagation. We prove numerically that this decay is only a transient 
feature corresponding to a reshaping of the input, and that a stable steady-state do exist. 
In other words we prove that the localized pulse of Ref. 46|] is really an ultrashort light 
bullet. This is the third result of our paper. 
Section IV presents our conclusions. 

II. DERIVATION OF THE TWO-DIMENSIONAL SINE-GORDON EQUATION 



We consider a set of two-level atoms with the Hamiltonian 

/ 



(1) 



v 



Ua 

Ub 

where f2 = Wf, — > is the frequency of the transition. The evolution of the electric field 
E is described by the wave equation 

(9J + 9f ) E = ^^d^ {E + AnP) , (2) 

where P is the polarization density. The light propagation is coupled with the medium by 
means of a dipolar electric momentum 

■•=(;:) 

directed along the same direction x as the electric field, according to 

H = Ho- fiE, (4) 

and the polarization density P along the x-direction is 

P = 7VTr(p/i), (5) 

where is the volume density of atoms and p the density matrix. Since, as shown in Ref. 



23| . the relaxation terms containing the relaxation times for both the level populations 



and coherences can be neglected, due to the fact that the relaxation occurs very slowly 
with regard to optical oscilations. Thus the density-matrix evolution equation (Schrodinger 
equation) reduces to 

ihdtp=[H,p]. (6) 

Transparency implies that the characteristic frequency of the considered radiation (in 
the optical range) strongly differs from the resonance frequency Vt of the atoms. Here we 
consider the so-called short-wave approximation, by assuming that Uyj is much larger than 
VL. 

Hence Vt is small, which can be expressed as 

Vt = en, (7) 

e being a small parameter. It is convenient to rescale the time unit in the basic equations 
according to ([7]), and hence we set 

c = ec, H = eH, t = i/e, (8) 

and analogously for Ua, Wb, -f^o- The set of equations (I!])-® is invariant under this transform, 
except that Eq. (j4]) becomes 

H = Ho- ^E. (9) 

It is equivalent to state that Q is small in the original units or that the typical wave frequency 
Uu), with respect to the rescaled time i, becomes very large, and corresponds thus to a fast 
variable 

Notice that we could equivalently have written t = t — z/V, as a zero order variable in the 
original units. However, the physical interpretation of the _approximation and its link to the 



standard short-wave approximation as defined in 45|, |47|, |48[ would be much less clear. The 
delayed time r involves propagation at some speed V to be determined. It is assumed to 
vary rapidly in time according to the assumption u^, ^ This motivates the introduction 
of the variables 

C = z, v = ^- (11) 
The pulse shape described by the variable r evolves more slowly in time, the corresponding 
scale being that of variable (. The transverse spatial variable y has an intermediate scale, 
comparable to what is usually considered in long- wave approximations j^ ]. 
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Next we use the reductive perturbation method as developed in Ref. 45|. To this aim 
we expand the electric field E as power series of a small parameter e: 

E = Eo + eEi+ e^E2 + e^E-s + .... (12) 

The polarization density P and the density matrix p are expanded in the same way. The 
expansion [Eqs. (fTT!) and (fT2|) ] is then reported into the basic equations ([I])-©, and solved 
order by order. The computation follows the same steps as the (l+l)-dimensional case, see 



Ref. 



23|. 



At order 1/e, the Schrodinger equation ([6]) reduces to 

ihdrPo = -Eq [p, po] . (13) 
We label the components of any Hermitian matrix u as 

u = [ . (14) 

\ U*t Ub J 

We assume that the coherences are zero long before the pulse, i.e., as z — > +oo or t — y — oo 

lim p,-4 = (15) 

T > — OO 

for any j > 0. Then the off-diagonal terms of Eq. f|T3l) allow to express the coherence po,t as 

Po,t = J[ J Eowodr', (16) 
where wq = po,6 — Po,a is the population inversion. We do not assume a pumping, hence 

lim Wo = wth, (17) 

r ^ — oo 

in which — 1 < wth < is the population inversion at thermodynamic equilibrium. If all 
atoms are initially in the fundamental state, then Wth = —1- The higher order terms of 
the population inversion [wj = pj^^ ~ Pj,a with j > 1) vanish at infinity. Using the above 
boundary conditions and definitions, we get from the diagonal terms of (IT^ the evolution 
equation for wq, as 

drWo = ^^^Eo I EoWodr'. (18) 
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Since the polarization density at this order is Pq = ^ (Po,*/^* + c.c.) = 0, where c.c. denotes 
the complex conjugate, we see from the wave equation ([2]) at leading order l/e^ that the 
velocity V should be = c up to this order. Eq. ^ at order e° is 

ihdrpi = Ho, po - Eq [/i, pi] - El [/i, po] • (19) 

The off-diagonal terms yield 

pi,t = iQ Po,tdT' + ^ / {EqWi + Eiwo) dr. (20) 

It allows to compute the leading term of the polarization density Pi = N {pi^tp* + c.c), as 



-2Nn \p 



2 



Pi = / / EoWodr"dT'. (21) 



oo ^ — oo 



Reporting Pi into the wave equation at order 1/e gives the evolution equation for the 
electric field 

which together with Eq. (fTSl) yields the sought nonlinear system of equations for Eq and Wq. 
Eqs. ffTSjl and fl22|) can be put in the normalized form 

BzT = AB + Byy, (23) 

At = -BC, (24) 

Ct = AB, (25) 

by setting 

y-f, ^=f. r = A^^, B = (26) 
with Lr a reference length in the micrometer 

and 

he , , 



The boundary conditions, taking into account Eq. flTTl) are then 



lim C = 0, lim A = — . (29) 

T — s>-oo T — s>-oo 
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The coupled system of partial differential equations (I23|) - (l25|) is a two-dimensional gener- 
alization of the sine-Gordon (sG) equation. It is worth to compare it to the set of equations 
derived by the same reductive perturbation method in the case of electromagnetic wave 
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-|44l|: however a careful analysis shows that system 



propagation in ferromagnetic media 

cannot be reduced to any of these equations. 
Next it is seen from Eqs. ( 124|) -( 125|) that Ot {A^ + C^) = 0. This allows us to introduce a 
function ip = ip^Z.T) as 

A = Ucosip, C = Usm.%lj, 



where U depends on Z only. From the boundary conditions we see that 



lim lb 

T — >-oo 







(30) 



(31) 



and thus U = Wth/wr, and is independent of Z (except if inhomogeneous pumping is present, 
but this situation is excluded here, see Eq. (refclt). Reporting Eq. ( l30l) into (HM and ( 125|) 
shows that ipx = B. Then, after integration using (IM]) . Eq. fl2^ reduces to 



ZT 



U sinip + ip- 



YYi 



(32) 



which is known as the two-dimensional sG (2D-sG) equation. Notice that, to the best of 
our knowledge, the 2D sG equation was never derived before by means of the short-wave 
formalism in any physical context. 

III. NUMERICAL RESOLUTION OF 2D SG 



The 2D sG equation is solved in the form fl23l)- (l25|) . using a simplified version of the 
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numerical scheme given in Refs 
plane wave with temporal shape as 

(T-Ti-To)2 



-|44j. Starting from an initial data in the form of a FCP 



B = /3exp 



Wn 



COS 



27r 



(T-Ti) + 7r/2 



in which 



Ti = 0.1 exp 



y2 



W 



Y 



(33) 



(34) 



yields a transverse perturbation. The parameters effectively used in computation were wt = 
1, wy = 2, To = —0.2, 6 = 1.5, f3 = 8, and U = —10. The numerical resolution clearly 
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Y Y 

FIG. 1: (Color online) Evolution of a perturbed input FCP plane wave into two dimensional FCP 
solitons. a) Input (Z = 0), b) Z = 52.2. 




FIG. 2: (color online) Two stages of the oscillations of the two dimensional FCP formed in the 
computation shown in Fig. [TJ a) Z = 51, h) Z = 52.2. 

shows that, during propagation, we get the formation of locahzed sohtons (see Fig. [1]). The 
two-dimensional FCP solitons are oscillating structures, which are localized in both space 
and time, see Fig. [2l 

The two-dimensional FCP solitons are roughly fitted by an expression of the form: 

B = Pexp[ J - — sin [uj{T - Ti)] , (35) 

in which the coefficients are chosen to fit the numerical data mentioned above. For the 
above data, this gives the following values: wt = 0.7, wy = 0.75, u = 3.5, (3 = —18, and 
Ti = 0.04. We see that the pulse stably propagates, being neither affected by dispersion nor 
by diffraction; after a transitory stage in which the pulse radiates energy, and its amplitude 
decreases, stabilization is reached eventually. The pulse energy is quite difficult to evaluate 
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due to the oscillations of the pulse. A pretty good evaluation of its evolution is found 
by considering the maximal value of the optical spectrum an its evolution versus Z . The 
corresponding plot is too noisy to be drawn here, but proves that stabilization occurs around 
Z = 100 for the particular data we used. The nonlinear inverse velocity varies slowly from 
0.65 to 1.08, and then stabilizes. The spectral width is also affected by the 'reshaping' and 
finally stabilizes, as shown in Fig. |H The spatial width of pulse is not modified at all, as 
can be seen in Fig. [5l Variations of the carrier-envelope phase are also observed, however 
they cannot be evaluated for technical reasons. 

In fact, the light bullets in the 2D sG model have already been mentioned in Ref. 46| . 



and even interactions have been studied in Ref. [491]. In 46|, the 2D sG equation was 
also derived from the Maxwell-Bloch equations, but the derivation was performed from a 
reduced form of the Maxwell-Bloch equations, and the physical assumptions were not so 
cleary given. We thus proved in the present work that the short wave approximation of the 
Maxwell-Bloch equations is indeed the 2D sG equation, and not one of the nonlinear system 
of equations 



Still in Ref. 



bund in the study of nonlinear electromagnetic waves in ferromagnets 
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43|. 



46|, a nonlinear Schrodinger-type equation with higher order nonlinear terms 



was derived from 2D sG, and it was shown that the latter saturates the nonlinearity and 
may arrest collapse. Light bullets were also simulated numerically in Ref. 46|], but on a 
rather short propagation distance; it was stated in that work that the light bullets loose 
energy and will become destroyed eventually. However, according to our computations on 



propagation distances 20 times larger than in Ref. 



461], in which stabilization of the energy 
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FIG. 4: (Color online) Evolution of the optical spectrum of the center of the pulse. 
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FIG. 5: (Color online) Evolution of the pulse width. The quantity maxr(i3)/ maxT^y(i?) is plotted 
versus Y for values of Z form to 317. 



occurs, the energy loss seems rather to be due to the reshaping of the input pulse, and the 
hght bullet is expected to have an infinite duration. Later, in Ref. 49|, the robustness of the 
light bullets during collisions was shown. The 2D sG equation was written in the form which 
allows propagation in both directions (in contrast, what we got here is the one-directional 
form of the 2D sG equation), and the inputs were either counterpropagating or colliding at 



a large angle 



49|. 



IV. CONCLUSIONS 



In summary, we have introduced a generic model beyond the slowly varying envelope ap- 
proximation of the nonlinear Schrodinger-type evolution equations, for describing the prop- 
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agation of (2+l)-dimensional spatiotemporal ultrashort optical solitons in two-level media. 
In the case of resonant nonlinearities, the short-wave approximation allow us to obtain a 
two-dimensional model generalizing the sine-Gordon equation. The use of the multiscale 
expansion up to the second-order in a certain small perturbation parameter allowed us to 
give a rigorous derivation of this equation, and more precisely to prove that the nonlinear 
partial differential equation obtained in this paper is different from the previously known 
two-dimensional generalizations of sine-Gordon equation, see, e.g., 4l|-|44|. Direct numerical 
simulations of the two-dimensional generalization of the sine-Gordon equation show that sta- 
ble (2+l)-dimensional ultrashort light bullets may form spontaneously from a transversely 
perturbed input plane wave, in the few-cycle regime. Thus the stability of light bullets has 
been demonstrated, in contrast with a previous reported result 

The present study can be generalized by taking into account both resonant and non- 
resonant optical nonlinearities. Moreover, the generalization to two transverse spatial di- 
mensions, in addition to time and longitudinal coordinates, in order to study the formation 
of (3+l)-dimensional few-optical-cycle light bullets (50| can also be envisaged. 
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